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Abstract 



The string-based Bern-Kosower rules provide an efficient way for obtaining parameter 
integral representations of the one- loop N - photon/gluon amplitudes involving a 
scalar, spinor or gluon loop, starting from a master formula and using a certain 
integration- by-parts ("IBP") procedure. Strassler observed that this algorithm also 
relates to gauge invariance, since it leads to the absorption of polarization vectors 
into field strength tensors. Here we present a systematic IBP algorithm that works 
for arbitrary N and leads to an integrand that is not only suitable for the application 
of the Bern-Kosower rules but also optimized with respect to gauge invariance. In 
the photon case this means manifest transversality at the integrand level, in the 
gluon case that a form factor decomposition of the amplitude into transversal and 
longitudinal parts is generated naturally by the IBP, without the necessity to consider 
the nonabelian Ward identities. Our algorithm is valid off-shell, and provides an 
extremely efficient way of calculating the one-loop one-particle-irreducible off-shell 
Green's functions ("vertices") in QCD. It can also be applied essentially unchanged 
to the one-loop gauge boson amplitudes in open string theory. In the abelian case, 
we study the systematics of the IBP also for the practically important case of the 
one-loop N - photon amplitudes in a constant field. 



1 Introduction 



It is by now well-known that techniques originally developed for the computation of amplitudes 
in string theory can be used also for the simplification of calculations in ordinary quantum field 
theory. Already in 1972 Gervais and Neveu observed that the field theory limit of string theory 
generates Feynman rules for Yang-Mills theory in a special gauge that has certain calculational 
advantages [lj. Actual calculations along these lines were done, however, only much later [21 
[3j HI [5] . A systematic investigation of the field theory limit at the tree and one- loop level was 
undertaken by Bern and Kosower, and led to the establishment of a new set of rules for the 
construction of the one-loop gluon amplitudes in QCD [Tl E| - These "Bern-Kosower rules" 
were used for a first calculation of the five - gluon amplitudes [9j. Their relation to the usual 
Feynman rules was clarified in |10j . 

Shortly afterwards, a simpler approach to the derivation of Bern-Kosower type formulas was 
initiated by Strassler |TTJ based on the representation of one-loop amplitudes in terms of first- 
quantized path integrals. For the case of QED representations of this type had been known for a 
long time [12J although they had rarely been considered as a tool for state-of-the-art calculations. 
Subsequently such representations were derived also for Yukawa and axial couplings P3J HH H5] , 
and generalized to higher loop orders [161 El OS 02] as wen as to the inclusion of constant 
external fields [2Ql ED E2] , gravitation ^ESESEHEIIEE] and finite temperature [251 1301 151] . 
For a review see 1321 . 



The central formula in the Bern-Kosower formalism is the following 'master formula': 
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(1.1) 



As it stands, this formula represents the one- loop N - photon amplitude in scalar QED, with 
photon momenta k^ and polarisation vectors £j. m denotes the mass, e the charge and T the 
total proper time of the scalar loop particle 0. 

Each of the integrals f dri represents one photon leg moving around the loop. The integrand is 
written in terms of the 'bosonic' worldline Green's function G B and its derivatives, 



Gb{ti,t 2 ) = \n-T 2 \-- 



T 

(n - r 2 ) 



G b (ti,t 2 ) = sign(ri - r 2 ) - 2 
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T 

(1.2) 



1 We work in the Euclidean throughout. With our conventions a Wick rotation kf — > —ik®, T — > is yields the 
N - photon amplitude in the conventions of |33| . 
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Dots generally denote a derivative acting on the first variable, Gb(t~i, T2) = ;^-Gb(ti, T2), and we 
abbreviate Gs%j = Gb(j% — T j) etc. In deriving the master formula a formal exponentiation has 
been used that needs to be undone by expanding out the exponential in (jl.ip in the polarization 
vectors, and keeping only the terms linear in each of them (in this paper our interest will be 
mostly in the off-shell case, but it will still be useful to keep the polarization vectors as book- 
keeping devices). 

We will not dwell here on the derivation of this formula, which can be obtained either from 
the infinite string tension limit of string theory [HI [3 El EI] or using the worldline path integral 
formalism [32] . Its role in the Bern-Kosower formalism is central, since it provides the input 
for the Bern-Kosower rules, which allow one to obtain, from the scalar QED integrand, the 
corresponding integrand for the photon amplitudes in fermion QED, as well as for the (on-shell) 
N - gluon amplitude in QCD. However, those rules do not apply to the master formula as it 
stands. Writing out the exponential in eq. (|l.ip one obtains an integrand 



exp<j • \ I multi-linear = (~i) N Pn{G Bij ,G Bij) exp 



1 X - 



2 ^ 

M =1 



TV 

(1.3) 



G Bijki • kj 



with a certain polynomial P/v depending on the various G Bij, G Bij and on the kinematic in- 
variants. The application of the Bern - Kosower rules requires one to now remove all second 
derivatives Gsij appearing in Pjy by suitable integrations by parts in the variables Tj. 

That this removal of all Cr^'s is possible for any N was shown in appendix B of [Tj- The 
new integrand is written in terms of the Gsij's and Gsij's alone, and serves as the input 
for the Bern - Kosower rules. Those allow one to classify the various contributions to the N 
- photon/gluon amplitude in terms of </> 3 - diagrams, and moreover lead to simple relations 
between the integrands for the scalar, spinor and gluon loop cases. A complete formulation of 
the Bern-Kosower rules is lengthy, and we refer the reader to [U [341 [32l [35] . For our present 
purposes, the most relevant part of the rules is that, up to a global factor of —2 correcting for 
the differences in degrees of freedom and statistics, the integrand for the spinor loop case can 
be obtained from the one for the scalar loop simply by replacing every closed cycle of Go's 
appearing in Qn according to the "replacement rule" 



G Bi\i^G Bi2i-i ' ' ' GBi„ii ^ G Bi^G Bi2i$ ' ' ' Gsi n ii G p{ x i 2 G pi 2 %^ • • • G ' Fi n ii 

(1.4) 

where Gfi2 = sign(ri — T2) denotes the 'fermionic' worldline Green's function. Note that an 
expression is considered a cycle already if it can be put into cycle form using the antisymmetry 
of Gb (e.g. GbyiGbyi = —Gbi2Gb2i)- A similar "cycle replacement rule" holds for the gluon 
loop case. Of course, in the nonabelian case there will also be many other modifications. 

As was discussed already in [7], the IBP procedure is generally ambiguous. However, this does 
not constitute an impediment to the application of the Bern-Kosower rules, whose application 
requires only that all second derivatives Gsij have been removed. For this reason, in the applica- 
tion to the computation of gluon amplitudes in [HI [9] the partial integration had been performed 
in an essentially arbitrary way. 

A closer look at the IBP was taken by Strassler in [36] . who noted that this procedure bears 
an interesting relation to gauge invariance. For each photon leg, define the corresponding field 
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strength tensor, 



Ff = m-sW (1.5) 

Remove all Gsij 's and combine all terms contributing to a given 'r - cycle' G bi^G Bi 2 iz • • • GBi n h ■ 
Then the sum of their Lorentz factors can be written as a 'Lorentz cycle' Z n (i\ii ■ ■ ■ i n ), defined 
by 



Z 2 (ij) 
Z n (hi 2 ...i n ) 

Thus Z n generalizes the familiar transversal projector. However, in [36] no systematic way was 
found to perform the partial integrations at arbitrary N, and also the absorption of polarisation 
vectors into field strength tensors (a process to be called "covariantization" in the following) 
worked only partially; after the IBP and the sorting of the resulting integrand in terms of "cycle 
content" some terms are just cycles or products of cycles, but, starting from the three-point 
case, there are also terms with left-overs, called "tails" in [36], and the polarisation vectors in 
them were not absorbed yet into field strength tensors. 

The IBP procedure was further studied in [37], where a definite partial integration prescription 
was given which works for any N, preserves the full permutation symmetry and is suitable for 
computerization. This algorithm is completely satisfactory as far as the application of the Bern- 
Kosower rules is concerned. However, it is obviously an interesting question whether some IBP 
algorithm exists which leads to an integrand where all polarisation vectors would be contained 
in field strength tensors, thus making gauge invariance, i.e. transversality, manifest at the 
integrand level. It is the purpose of the present work to present such an algorithm 0. 

In detail, we will do the following: In chapters [2] and [3] we find a surprisingly simple way of using 
IBPs to covariantize the Bern-Kosower master formula itself; the resulting representation will be 
called the R-representation. In|3]we summarize the "symmetric IBP procedure" of [37] . leading 
to the Q - and Q' - representations. The next two sections [5] and [6] define our new algorithm, 
which combines elements of both the R - and the Q' - representation, and results in what we 
will call the S-representation of the N photon/gluon amplitudes. As an aside, in section [7] we 
present a further improvement of the "two-tail" which leads to a particularly compact integrand 
at the four-point level (but does not seem to generalize to the N - point case). In section [8] we 
shortly comment on a direct treatment of the spinor QED case in the worldline super formalism. 
Section [9] is devoted to our main application, which is the calculation of the one-loop off-shell 
one-particle- irreducible N - gluon amplitudes (or "N- vertices" ) . While in the abelian case there 
are never any boundary terms in our IBPs, since all integrations run over the full loop and 
the integrand is written in terms of worldline Green's functions with the appropriate boundary 
conditions, in the nonabelian case the color ordering of the gluon legs leads to the restriction 
of the multiple parameter integrals to ordered sectors, and to the appearance of such boundary 

2 An IBP algorithm for the general worldline integrand was also developed in [38], but for the unrelated purpose 
of tensor reduction. 
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(1.6) 
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terms [11] 136]. Those generally can be combined into color commutators and, in x-space, would 
in principle allow one to achieve a complete nonabelian extension of the covariantization, namely 
to rewrite the final integrand in terms of full nonabelian field strength tensors, and to complete 
all derivatives to covariant ones. This is not possible in momentum space, but here instead 
the IBP procedure generates a natural form factor decomposition of the iV-vertices, where the 
bulk terms are manifestly transversal and all non-transversality has been pushed into boundary 
terms. Finding such a decomposition by standard methods usually involves a tedious analysis 
of the nonabelian Ward identities, and so far has been completed only for the three-point case 
[39]. 



In section [TU] we return to the abelian case, but now with a constant external field added. Here 
the replacement rule (jl.4p applies as well [21] [22] , however the IBP procedure is complicated 
by the fact that the worldline Green's functions become nontrivial Lorentz matrices. We will 
point out here the necessary modifications, and also find a way of using the IBP procedure to 
effectively eliminate the nonvanishing coincidence limits of the generalized Green's functions Qb, 
Qf (see (|10.2p . (|10.5p below) which otherwise would lead to a proliferation of terms as compared 
to the vacuum case. 

In the conclusions section we summarize the properties of our new IBP algorithm, and discuss 
various applications, some of which have already been published or are actually in progress. We 
shortly comment on possible generalizations to gravity and string theory. 



2 The P - representation 

We will call "P-representation" the integrand obtained directly from the expansion of the Bern- 
Kosower master formula, 



T scal [ki,e i; ...;k N ,e N ] = (-ief (2tt) d 5(J2 h) J ^AnfT)-^ e~ m " T 

N r T ■ r N i i 

x II / dT i p N(GBij,G B ij) exp< ^2 2 G Bijh ■ kj \ 
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(2.1) 



Explicitly, the polynomial P/y is given by 

Pn = Gsih ei-k^ Gb2i 2 £ 2 • h 2 ■ ■ ■ Gbni n £ n ■ h N 

N ______ 

— ^2 GBa,b £ a-£bG Blh^l' &ii ' ' ' GBaiSa-k{ a ■ ■ ■ G ' Bbi b ^b'ki b • • • GBNi N £ N'ki N 

a, 6=1 
a<b 

N 

+ ^ [G B ab £ a ' £ bGscd £ c m £ d + GBac £ a ' ^cGsbd^b ' £ d + Gsad £ a ' £ dGsbc £ b ' £ c) 

a,b,c,d—l 
a<b<c<d 

XGBHiSl-kii ' ' ' G Bai a ^a' h a ' ' ' GBbi b ^b'ki b ' ' ' GBci c ^c'ki c ' ' ' GBdi d £ d'h d ' ' ' GbN{ n £ N ' &i N 

(2.2) 
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Here and in the following the dummy indices 11,12, ■ ■ ■ should be summed over from 1 to N, 
and a 'hat' denotes omission. Note that all terms in Pjv are obtained from the first one by a 
simultaneous replacement of pairs of GBri r £ r • h r GBsi s £s • h 3 by —Gsrs^r • £s, which has to be 
done in all possible ways. Note also that Gbu = by antisymmetry. 

These P-representation integrals are still directly related to the ones arising in a standard Feyn- 
man or Schwinger parameter calculation of the N photon amplitude [IP} [TT] . The exponential 
factor will, after a multiple rescaling and performance of the global T - integration, turn into 
the standard one-loop N - point Feynman denominator polynomial. The 5 - function contained 
in Gsij will bring together the photons i and j, corresponding to a quartic vertex, and the 
contributions of such terms match the ones from the seagull vertex of scalar QED. 

3 The R - representation 

Before coming to the "old" IBP procedure of 137] and its intended improvement, it will be 
useful to solve a simpler problem, namely how to covariantize the Bern-Kosower master formula 
itself. In the following we will often abbreviate 

e ^T, G Bijk i -k j _ e (-) / 3> -j\ 
Consider first the case of N = 2, where 

Pi = G B i2£l ■ k 2 G B 2i£2 ■ h - G,Bl2£l ■ £2 (3.2) 
We choose two vectors r\ , T2 that fulfill 



n ■ h ^ (3.3) 

but are arbitrary otherwise. Adding to the integrand P2 e*-'-* the following sum of total derivative 
terms, 



ri ' £l di (g B 2iE2 ■ h e«) - r -^d 2 (G m2 ei ■ k 2 e«) + ^f^^ 2 e« (3.4) 



ri • ki V / r 2 ■ k 2 V J r\-k\r 2 -k 2 

{di = -Mr) the total result is a change of P2 into R2, 



„ . a ri • F\ ■ k 2 ■ r 2 -F 2 -ki ~ r l -F l -F 2 - r 2 

til •= ^B12 ; <~rB2l ; h <-^_Bl2 ; ; — (0.0) 

n ■ ki r 2 - k 2 ri ■ kir 2 • k 2 

Thus we have managed to absorb the polarization vectors into field strength tensors. And this 
procedure can be immediately generalized to the Appoint case: let us abbreviate 
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T r (l) 



r i ' Fj ■ kj 

G B i 



W r (ij) := G Bt 



r i " kj 



f"i ' Fi • Fj ■ Tj 



(3.7) 
(3.8) 



and choose vectors j-j., . . . , rjy fulfilling (j3.3[) . Then it is a matter of simple combinatorics to 
verify that 
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A (> ' + [n^-Pa^AJ-l 

a=l 
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where the operator A a is defined as follows: each term in P/v either involves the index a 
in a second derivative factor Gb, or it carries a factor of Gs a i a e a ■ ki a . In the former case the 
term will be annihilated by A a , in the latter case the action of A a is to replace the factor of 
G B ai a e a ■ h a by 1. 

We can then reexponentiate the new integrand, and arrive at the following covariantized version 
of the Bern-Kosower master formula (jl.ip : 



N rT 
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Thus we have achieved manifest gauge invariance at the integrand level, with a large freedom 
of choosing the vectors r±, . . . ,rjy. We will call this the "R-representation" of the ./V - photon 
amplitudes. Note that it reduces to the original master formula (jl.ip if r i ■ e% = for all i. 



4 The Q and Q' - representations 

Next, we review the IBP procedure motivated by the Bern-Kosower rules, whose primary purpose 
is to get rid of all second derivatives Gb [2 GEE 157] . 

For the two-point case P2 has been written down already in (|3.2 j) . After an IBP of the second 
term in either t\ or t 2 , and using Gbyi = —Gb21, it turns into 



Q2 = Gbi 2 Gb21 (ei • k 2 e 2 ■ k x - s x ■ e 2 h ■ k 2 ) = G B i 2 G B 2iZ 2 {12) (4.1) 



Proceeding to the three-point case, here (|2.2p becomes 



P3 — Gbh£\ • hGB2j£2 ■ kjGB3k^3 m kk 

- GbuSi ■ e 2 G B 3i£3 ■ h + (1 -> 2 -)• 3) + (1 -> 3 -> 2) 



(4.2) 



In this three-point case it is still possible to remove all G^'s in a single step. To remove, e. g., 
the term involving GbwGbzi m the second term of P3, we can add the total derivative term 



#2 ( G B 12Sl ■ £2Gb31^3 • ^ie ( ) ) (4.3) 



This term together with five similar ones removes all the Gb's. Decomposing the new integrand 
according to its "cycle content", P3 gets replaced by Q3 = Q3 + Q|, where 



Q3 = Gb12Gb23Gb3iZ3(123) , 

Ql = Gbi2Gb2iZ 2 (12)T(3) + Gbi3G B 3iZ2(13)T(2) + G B23 G B 32^2(23)r(l) 

(4.4) 

Note that Q3 contains a cycle of length three and Q| a cycle of length two, as indicated by the 
upper indices, and that each r-cycle appears together with the corresponding "Lorentz-cycle" . 
This motivates the further definition of a "bicycle" as the product of the two: 



G(i\i2 ■ ■ - in) '■= GBi 1 i 2 GBi 2 i 3 ' ' ' GBi n i x Z n (i\l2 ■ ■ - in) (4.5) 

But the terms of Q% have, apart from the cycle, also a "one-tail" , defined by 



T{a) := G Bm £a ■ h (4.6) 

This tail still has a polarisation vector that is not absorbed into a field strength tensor. It is 
easy to see that, nonetheless, each term in Q3 is individually gauge invariant; if one replaces in, 
e.g., the term 

G B 12GB2lZ 2 (12)GB3kS3-he^ 

£3 by &3, then it becomes proportional to 



5 3 (GBi2G B2 iZ 2 (12)e») 

However, our aim here is to make gauge invariance manifest even at the integrand level. Now in 
the three-point case there are already various chains of IBP that can be used to remove all the 
Gb's, but if one assumes that the corresponding total derivative terms are added with constant 



S 



coefficients (i.e., they involve no dependences on momentum or polarization other than the ones 
already present in the term which one wishes to modify) , then it is easy to convince oneself that 
they all lead to the same Q 3 of (|4.4jl . Thus we have to look for a more general type of IBP. We 
will now essentially apply the procedure of the previous section to the tails. Consider again the 
first term in Q\ above, eq. (|4.4j) . Choose a momentum vector r 3 such that r 3 • k 3 7^ 0, and add 
the total derivative 



Z 2 {12)d 3 [G Bl2 G B21 e^) (4.7) 



_ 7-3 • £3 
?"3 ■ h 

The addition of this term to the first term in Q3, and of similar terms to the second and third 
one, transforms Q 3 into 



?2 ._ n n v mo\/» r 3 ■ F 3 ■ k k ■ ■ ■ r 2 ■ F 2 ■ kj 



R 3 :— G b\ 2 G b 2 iZ 2 (12)G Bzk ; \- GBizGBz\Z 2 {l2>)GB2j 



+ G B23 G B32 Z 2 (2?>)G B i 



r 3 -k 3 r 2 ■ k 2 

i 



n ■ ki 



(4.8) 



Thus we have completed the covariantization of the integrand. 



In the abelian case the 3-point amplitude must, of course, vanish, which we can see by noting 
that the integrand is odd under the orientation-reversing transformation of variables Tj = T — r-, 
i = 1,2,3. 

Proceeding to the four-point case, here even using only total derivative terms with constant 
coefficients (in the above sense) there are already many ways to remove the G^'s by IBP, 
with a large ambiguity for the final integrand, and it is not obvious how one should proceed. 
But certainly one would like to preserve the manifest permutation symmetry of the iV-photon 
amplitudes, and in [37] this was used as a guiding principle to develop the following algorithm: 



1. In every step, partially integrate away all the second derivative factors Gsi/s appearing 
in the term under inspection simultaneously. This is possible since different Gb^'s never 
share variables. 

2. In the first step, for every factor of G B ij present use both Tj and Tj for the IBP, and take 
the mean of the results. 

3. At every following step, any Gsij appearing must have been created in the previous step. 
Therefore either both variables Tj and Tj were used in the previous step, or just one of 
them. If both were used, then both should be used again in the actual IBP step, and the 
mean of the results be taken. If only one of the variables was used in the previous step, 
then the other variable should be used in the actual step. 



In the four-point case, applying this algorithm to P4 and decomposing the resulting integrand 
according to cycle content, leads to the following version of the numerator polynomial Q4 |37| : 
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Qa — q\ + q\ + q\ — Ql 2 

Q\ = G(1234) + G(1243) + G(1324) 

Q\ = G(123)T(4) + G(234)T(1) + G(341)T(2) + G(412)T(3) 

Q\ = G(12)T(34) + G(13)T(24) + G(14)T(23) + G(23)T(14) + G(24)T(13) + G(34)T(12) 
Qf = G(12)G(34) + G(13)G(24) + G(14)G(23) 

(4.9) 

Here we have now further introduced the two-tail, 



T^j) '■= G B ir£i ■ k r G B js£j ■ k s + -G B ij£j ■ £j G B jrh ■ K - G B j r kj ■ k r 1 (4.10) 

r,s ^ ~ ' 

Thus the final representation of the (still off-shell) four-photon amplitude in scalar QED becomes 



4 pOO 7^7"! 

r sca i[k 1 ,e 1 ;...;h,e4 = -^(2i)Wy —T t ( ' 

(4tt)^ ^ Jo T 

x J dr x --- dT A Qi{G B ij) exp j ^ ^ G B ijh ■ kj | 



(4.11) 



We shortly summarize the advantages of this representation compared to a standard Feyn- 
man/Schwinger parameter integral representation (see |32j for details): 

First, the rhs of (|4.1ip represents already the complete amplitude, with no need to add "crossed" 
terms. The summation over "crossed" diagrams which would have to be done in a standard field 
theory calculation here is implicit in the integration over the various ordered sectors. 

Second, the IBP procedure has homogenized the integrand; every term in Qn has N factors 
of G B ij and N factors of external momentum. In the four-point case this has the additional 
advantage of making the UV finiteness of the photon-photon scattering amplitude manifest 
before integration. While the original numerator P4 contains terms involving products of two 
G B ij's which lead to spurious divergences in the T - integration, after the IBP the integrand is 
finite term by term. 

Third, It allows one to obtain the corresponding spinor QED amplitude by the application of 
the replacement rule ()1.4|) . In applying the rules it must be observed, though, that the form of 
the integrand given aboven still contains, apart from the explicit cycle factors, additional cycles 
from the tail factors for certain values of the dummy indices. In the four - point case this occurs 
for Q\ only: The two - tails contained in Q\ as given in (|4.9p above each contain a two-cycle, 
since the content of the braces on the rhs of (14. 101) for r = j, s = i turns into 



G B ijG Bji (ei ■ kjEj ■ ki- Si - Ejh ■ kA = G(ij) (4.12) 



10 



For the application of the "replacement rules" it is therefore convenient to decompose Q4 in a 
slightly different way [37]. Namely, note that 



Ql = Qi + 2Qf (4.13) 

where Ql is obtained from Q\ by eliminating the term with r = j, s = i from the sum over 

dummy indices. With this definition, and setting Q 4 = Q\ for the remaining components, we 
can write 



Q 4 = Qt + Qi + Qi + Qt (4.14) 

In this form all cycle factors are explicit, and moreover all the coefficients in the decomposition 
turn out to be unity. Generally, we will denote by T'(i\ . . . i n ) a tail whose cycles have been 
removed by the appropriate restrictions on the multiple dummy sums appearing in it. 

Thus the four-photon amplitude for the spinor loop case can now be obtained from the scalar 
loop formula (14. lip simply by multiplying with a global factor of —2 from spin and statistics, 
and by replacing, simultaneously, each bicycle G{i\ . . . i n ) by the corresponding "super-bicycle" 



Gs(h ■ ■ - in) '■— {G Bhi 2 G Bi 2 ii ' ' ' G Bi n i x ~ G Fhi 2 G Fi 2 i 3 ' ' ' Gpi n i x )Z n (ix ■ ■ ■ i n ) (4.15) 

(the notation refers to the worldline supersymmetry underlying the replacement rule f)1.4|) . see 
|32j). This "symmetric partial integration" procedure has been worked out explicitly for up to 
the six-photon case; see [32j [37] for the explicit formulas. 



5 The QR representation 

So far we have established two seemingly unrelated IBP procedures, the first one leading to the 
manifestly gauge invariant R-representation, the second one to the Q' - representation that is 
suitable for the application of the Bern-Kosower rules, but manifestly gauge invariant only in the 
cycle factors, not in the tails. We will now combine the two IBP strategies, using the following 
three simple observations: 



First, it had been noted in the appendix C of [32J that the Q - representation is recursive, in the 
following sense: Each term in the cycle decomposition of Qn contains at least one cycle [32], so 
that any tail appearing in the ./V - point amplitude has at most N — 2 arguments. And a tail of 
length, say, M, is related to the (undecomposed) lower-order Qm simply by writing Qm in the 
tail variables, and then extending the range of all dummy variables occurring in it to run over 
the full set of indices 1, . . . , N. For example, writing out (|4,10p for the two-tail T(12) in the last 
term of Q\ in (14. 9p gives 



4 4 
T(12) = ^2 G B \r£i ■ k r G B 2s£2 ■ k s + -G B i2S\ ■ e 2 ^ [GW&i • K - G B 2rk% ■ k r (5.1) 

r,s=l ~ r=l 
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which can also be obtained by writing Q 2 , defined in (|4.ip . as 



Q 2 (12) = G B i 2 ei • k 2 G B21 e 2 • h + -G B12 si ■ £2 



G B \ 2 ki ■ k 2 - G B2 \k 2 ■ ki 



(5.2) 



and introducing appropriate dummy index summations. That this property holds in general can 
be easily seen by considering those terms in the cycle factors of the decomposition of Q m that 
do not contain factors of £{ ■ £j, and thus can have involved IBPs only in the tail and not in the 
cycle variables; see the appendix C of [32] for more details. 

Second, for each N the symmetric IBP procedure defines a unique Qm and thus a total derivative 
term 



S N e^ =Q N eU-P N eU 



(5.3) 



Consider now an arbitrary term in the cycle decomposition of Qn- It will have the form 
C(i\ . . . iL)T{ji . . . jm), where M < N—2, C(i\ . . . il) is a bicycle or product of bicycles involving 
the variables t« 1 , . . . ,n L , and T(ji . . . jm) is the unique (in the symmetric IBP scheme) tail of M 
variables, written in the remaining variables Tj x , . . . , tj m . Consider C(i\ . . . i^SmUi ■ ■ ■ 3m) 
where possible dummy variable summations in Sm are extended to run over the full range of 
variables ti,...,tjv as above. This is still a total derivative term (involving only derivatives in 
the tail variables), and the above simple relation between the M - tail and Qm implies, that 



C(h . . . i L )T(ji . . . j M ) e (0 - C(h . . . % l )SmUi ■ ■ ■ 3M) e° = C(h . . . i L )T p {h . . . j M ) e ( ' } (5.4) 

with a new version T p of the M - tail which relates to Pm in the same way as the standard tail 
T to Qm, by an extension of the dummy index sums. 

Continuing with our example above, here we have 



and (15.41) becomes 



(Q 2 - J P 2 )e Gsi2fclfc2 



12^1 ' ^2 



Gswki ■ k 2 - G B 2\k2 ■ h 



l(d 1 -d 2 )(G B12 e 1 -e 2 e G ^ k ^ 



+ G B12 e 1 -e 2 }e G ™ k ^ 



(5.5) 



G(34)T(12) e« -\{di~ d 2 ) (G(34)G B12£l • e 2 e«) = G(34)T P (12) e« (5.6) 



where T(12) was given in (|5.ip and T p (12) is given by 



4 

T p (12) = GsirSi ■ k r G B2s e 2 ■ k s - GbuE! ■ e 2 (5.7) 
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Finally, the new tail T p (- ■ ■ ) can now be covariantized by a simple extension of (13. 9h : 



- p ia 5 ia A j j(c(-)T p (j 1 • ..j M ) e«) = C(-)T p (gW« • fc ia -»• T r (a), -<Wa • £6 -> W r (a6)) 

a=l 



In our two-tail example, this lead to the following result: 



T r (12) = G B1 /-±^G B2 /-^^ + G B12 ^4-^ (5.9) 
where the sums over r and s run from 1 to 4. This should be compared with R 2 of (|3.5|) . 



6 The S representation 



Thus in the QR - representation we have the usual decomposition into cycle and tail factors, with 
the tails already covariantized, and in a form that generalizes the (lower order) R - representation 
by an extension of the dummy index sums to run over all N variables, including those belonging 
to the cycle factors of the term under consideration. To finish our quest for an integrand that 
would be both covariant and suitable for the application of the Bern-Kosower rules, two more 
steps are needed: First, we need to remove the remaining Go's; this can be done by reapplying 
the symmetric IBP procedure of chapter 01 without any modifications. And finally all cycles still 
contained in the tails have to be separated out. We will call this final, in some sense optimized 
result for the integrand of the N photon/gluon amplitudes, the "S - representation", and denote 
the corresponding tails by T s '(-). 

Continuing with our example of the two-tail in Q|, the first step transforms T r (12) of (|5.9p into 



21(12) - £ G Bl /-^^d ms ^-^ 
^ n • k\ r 2 ■ k 2 

r,s=l 

4 



t:Gbi2 (y~] Gsirh ■kr-'Y] G B 2sk 2 ■ k s ) Tl F ] r2 
2 \f— f \ I n • kik 2 ■ r 2 

(6.1) 



r=l s=l 



This form of the two-tail, like the original two-tail T(12) of (|4.10p . still contains a cycle - the 
terms on the rhs with r = 2 and s = 1 combine to form a G(12), as in (|4.12p . Eliminating these 
terms from the tail one arrives at the final form, 
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r s '(12) . Yl G B1 /-^^G m /-^-^ 
^ T\-k\ r 2 • k 2 

~G\B12 (y~) G B lrh -kr-Y] G B 2sh ■ k s ) r l F l F 2 r 2 
; V- 2 -^ J T\-k\k 2 -T 2 



r,s = l 
r-,s#(2,l) 



r=l 
r^2 



s = l 



.2) 



For the one-tail there is no difference between T r (i), T a (i) and T s '(i), being all given by (|3.7p . 

We can now write down a covariantized version of the four photon amplitude, simply by taking 
over (|PJ) . (|4.1ip . and (|4.14|) and replacing all one-tails by ([3"7T|) and all two-tails by (|6.2p . 
Explicit formulas for higher point amplitudes will be given elsewhere. 



7 Alternative version of the two-tail 

For the two-tail, there is actually yet another form which is covariant, free of Gb's and at the 
same time more compact than f|6. 1 j) . Starting again with T p (12) of ()5.7[) we add the following 
total derivative term to T p (12) (omitting now the inert cycle factors C(-)) 



' tr(FiF 2 )ai9 2 e () + ' 



(A* • k 2 ) 2 



h ■ k 2 - 



e\ ■ £ 2 <9i<9 2 e () - e\ ■ k 2 e 2 ■ kjdi ( G B2j e () ] - e 2 ■ k\£\ ■ kid 2 \G B \ 



One obtains the new two-tail 



Iff (12) = G B nG B2 jki ■ H12 ■ kj 
where we have introduced the tensor 



(7.2) 



"12 — 



(F 1 F 2 yk 1 ■ k 2 - fc£fc|tr (FiF 2 



(7.3) 



Note that tri/i 2 = and Hj 2 = H 2 \. Note also that the term with i = 2,j = 1 in (|7,2p as 
before produces a G(12), since 



k 2 ■ H 12 ■ h = itr (FxF 2 ) - tr (F 1 F 2 ) = ~tr {F 1 F 2 ) = -Z 2 (12) 



(7.4) 



Thus Tjj(ij) can be used as well as T(ij) and T s (ij) in the construction of the four-point 
amplitudes, including the application of the replacement rule (|1.4|) and the simple sign change 
in passing from (|4.9p to (|4.14p . However, contrary to T s (ij) it appears that Tu(ij) has no natural 
generalization to the higher-point tails. 
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8 The case of spinor QED 



One of the main purposes of the IBP procedure is to trivialize the transition to the spinor QED 
case though the replacement rule (jl.4p . Still, it is interesting to note (and of possible practical 
relevance for the generalization to the case of open fermion lines, where the IBP is less attractive 
due to the existence of boundary terms) that in the worldline formalism there is also a more 
direct treatment of the spin ^ case using an approach based on explicit worldline supersymmetry 
\40\ \T7\ \32\ [38l |4"T] . It allows one to write down a master formula for N - photon scattering |40j 
which is formally analogous to the one for the scalar loop, eq. fll.ip : 



poo jrp D 

T spin [k 1 ,e 1 ;...;k N ,e N ] = -2(-ie) N (2tt) d 5( V h) / — (^Ty^e 

Jo 1 

N T , N 

x n / dn I de i eyi p\ ^2 

i=l ® "' i,j=l 



D_ 2 



m 2 T 



1 . . 1 

~Gijki ■ kj + iDiGijEi ■ kj + —DiDjGijEi ■ £j 



lin(£i,...,e N ) 



Here we have further introduced integrals over the Grassmann variables 9\, . . . ,9n, such that 
J dOiOi = 1, and the super derivative 



The two worldline Green's functions Gb,f now appear combined in the super Green's function 



6(ti,0i;75,0 2 ) = G b (t 1 ,t 2 ) + 6 1 6 2 G f (t 1 ,t 2 ) (8.3) 

Now also the polarization vectors £i,...,£jv are to be treated as Grassmann variables. The 
overall sign of the master formula refers to the standard ordering of the polarization vectors 
eis 2 ■ ■ ■ en- 

Starting with this master formula, all the manipulations which we have applied in the previous 
chapters to the scalar loop integrands can, mutatis mutandis, also be used in the spinor loop 
case starting from (|8.1|) . Here we will be satisfied with pointing out that the covariantized 
Bern-Kosower master formula (|3.10p generalizes to the spinor QED case as follows: 



pOO jrp 

T spiQ [k 1 ,e 1 ;...;k N ,e N ] = -2(-ie) N (2tt) d 5(Y / h) J jri^T) 

N rT „ , N 



i=l 



n / dr i I rf ^exp<j 



-Gijh ■ kj + iDidj °- -DiDjGij 3 3 



lin(Fi,...,F N ) 
(8.4) 



where now Fx,... , Fn have to be treated as Grassmann variables. 
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9 The nonabelian case 



As far as concerns the calculation of the on-shell QED photon amplitudes, or of the on-shell 
gluon S-matrix elements via the Bern-Kosower rules, the availability of a representation that 
is manifestly transversal at the integrand level is satisfying, but the significance of this fact for 
practical calculations is not obvious. To the contrary, it is easy to recognize the advantages of 
such a representation when it comes to the gluon amplitudes off-shell. Here the natural objects 
to consider in QCD are the "N- vertices" , that is the one-particle-irreducible N - point functions, 
and for applications of those it is often essential to decompose them into a basis of transversal 
and longitudinal tensor structures (see, e.g., [391 S21 S31 S3] ) - Such a "transversality-based form 
factor decomposition" in the present approach emerges essentially automatically in the IBP 
procedure through the appearance of field strength tensors. We have seen how this happens for 
the abelian case, but it is true also for the nonabelian case; here in principle one would like to 
see the full nonabelian field strength tensor emerging, 

F, v = F^T a = F® u + ig[AlT\ A C U T C ] (9.1) 

where by 



F%, ee {d,Al - d„A«)T a (9.2) 

we now denote its "abelian part"; and indeed Strassler demonstrated already for some simple 
cases how this happens 136] : When the external particles are gluons, the various ordered 
sectors of the integral J dr\ ■ ■ ■ dT n need to be considered separately, since they carry different 
color factors. Therefore boundary terms now arise in the IBP procedure, and the commutator 
terms are generated as differences of boundary terms between adjacent sectors that in the abelian 
case would cancel, but cannot do so any more in the presence of color since two of the color 
matrices appear in different orders. In an x-space calculation of the effective action, those 
commutator terms could then be combined with the "abelian" parts of the field strength tensor, 
but this is not possible in a momentum space calculation of the N - point function at fixed 
N, since any term in the nonabelian effective action after Fourier transformation contributes 
to amplitudes with various numbers of external particles; e.g., the term tr (D^FapD^F 01 ") will 
contribute to the N - point functions with N between two and six. Generally, each term in 
the nonabelian effective Lagrangian has a "core" term, which has a counterpart already in the 
abelian case (in the example this would be d ll F^ /3 d' J 'F 0a ^) and a number of "covariantizing" 
terms that all involve commutators, and belong to amplitudes with more legs than the core 
term. In this section, we will explain the essentials of how to calculate the scalar, spinor and 
gluon loop contributions to the one-loop N - gluon vertex. The details and a full recalculation 
of the three-gluon- vertex will be left to a separate paper |45j . 

Starting with the scalar loop case, here the master formula (jl.ip generalizes to the nonabelian 
case simply by supplying a global color factor, and keeping the gluons in a fixed order: 
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/■oo 

r^^[fci,ei;. = Hg) N ix(T a K..T a »)(2n) D i8(Y l ki) dT(^T)- D ' 2 e 

J o 

>>T /-ri /"Ov-2 f JX^ 

x / dri / dr 2 . . . dr N _i exp < >J 

Jo Jo Jo ■ • i 



71 '" r i /"tjv-2 ( rx . x 



-m 2 T 



lin(ei, 



( 



Here the T a are the generators of the gauge group in the representation of the loop scalar. This 
treatment of color correponds to a "color-ordered" representation (although not necessarily in 
the usual sense, where the T a would be in the fundamental representation, see, e.g., (46J 147]). 
Note that we have not only fixed the ordering of the gluons along the loop but also used the 
translation invariance in r to set r^r = 0. Summation over all (JV — 1)! inequivalent orderings of 
the iV vertex operators is implied. Also, T has been given an index "1PI" to indicate that the 
rhs gives only the one-particle-irreducible part of the N - gluon amplitude, not including the 
reducible part that now also exists, differently from the abelian case. Starting from (|9.3|) one 
can apply the IBP procedure leading from the P - representation to the S - representation as 
before, the only novelty being the boundary terms. Since for the bulk terms all the polarization 
vectors £j get absorbed into the transversal structures (jl.5p (which now, however, represent only 
the "abelian part" of the field strength tensor), in the final representation the non-transversal 
part of the N - vertex must be entirely in the boundary terms, given by lower-point integrals. 
For those one still has to choose the vectors r^, preferably in a way that is consistent with the 
cyclic invariance of the nonabelian amplitudes. In the three-point case a convenient cyclic choice 
is n = &2 — &3> r 2 = &3 — k%, r 3 = k\ — k2, and indeed it turns out [15] that, with this choice, 
the resulting form factor decomposition matches precisely the standard Ball-Chiu decomposition 
|39j of the three-gluon-vertex. 



Coming to the spinor loop case, here the only issue is whether the replacement rule (|1.4p can 
be applied also to all the boundary terms now arising in the IBP procedure. This is indeed the 
case, as we can see as follows: it suffices to show the corresponding statement for the effective 
action, rather than the momentum space Green's functions. Now, the effective Lagrangian can 
in principle be written as an infinite series of terms that are Lorentz scalars formed using any 
number of field strength tensors and covariant derivatives. As was already mentioned, each such 
term has a core term, whose calculation is not different from the abelian case for either the 
scalar or spinor loop, such that the replacement rule applies to it. All covariantizing terms of a 
core term must share its coefficient, and low-order calculations show that, as one would expect, 
the way this works is that they all involve the same parameter integral [36j [45| . And for the 
whole structure to continue to be gauge invariant for the spinor loop case it is necessary that 
the same replacement rule applies to all the covariantizing terms as well as for the core term. 

Similarly one can convince oneself that also the replacement rules that connect the scalar with 
the gluon loop cases [El [Ml [32], [Til EI] can be extended from the core terms to the ones involving 
boundary contributions |45| . An additional issue with the gluon loop contribution to the ./V - 
vertex is that one has to choose a gauge for the gluon propagator. The application of the gluonic 
replacement rules gives the N - vertex corresponding to the use of the background field method 
with Feynman gauge for the quantum part [21~] , which is also known to coincide with the 
result of the application of the pinch technique [48|, |4"9T [50] . This version of the gluonic vertex 
is also the one that leads to SUSY sum rules [5T] . 
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10 The constant external field case 



The Bern-Kosower master formula (jl.ip has the following straightforward extension to the 
N— photon amplitude in a constant external field \20 \ I21|, [22] : 



T scal [h,e 1 ;...;k N ,e N ] = (-ie) JV (2vr) D <5(^ h) / [4vrT] — e-^det^ 



o 



sin Z 



i N r T 

n ' dr < 



i=l 



' ^-Bij ' ' Gsij ' kj + — £j • G Bij ' £j 



llin(£i,....e N ) 

(10.1) 



Here we have introduced the matrix Z^ = eTF^ , Gsij = GBiji^j) denotes the "bosonic" 
worldline Green's function in the constant field background, and Gsi^GBij its first and second 
derivatives with respect to Tj. Explicitly, they are given by I22j 



Qbi2 
Qbyi 
Qb\2 



T ( Z 



2Z 2 Vsin-E 
/ / Z 



e -iZG B12 + iZ G m2 _ i 



Z V sin Z 



2S(n - r 2 ) 



%ZGb\2 



2 Z 
T sin Z 



-iZGbi2 



while the fermionic Green's function Gpij generalizes to 



(10.2) 



SF12 = Gf 



12" 



—iZGb\2 

cos Z 



(10.3) 



These expressions are to be understood as power series in the matrix Z. We note the symmetry 
properties 



GBji = GBij, GBji = -GBij, GBji = GBij, <3Fji = -GFij (10.4) 



and the coincidence limits 



Gb(t,t) = ^L(ZcotZ-l) 

% 

Gb(t,t) = icotZ - — 
Gf(t,t) = —it&nZ 



[10.5) 
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Note that those are independent of r. For the following, it will also be convenient to introduce 
the "subtracted" Green's function 

T ( e~ iZ( ^ B12 — cos Z ■ \ 

Qb(ti,t 2 ) := Gb(ti,t 2 ) -Gb(t,t) = — ^— ViG B \2 (10.6) 

LZ \ smZ / 

The corresponding spinor QED amplitude is obtained from this master formula in the same way 
as in the vacuum case, with the following minor modifications: First, note from (|10,2p . (|10.5p that 
the worldline Green's functions for the constant field are generally nontrivial Lorentz matrices. 
Thus the definition of a 'cycle' must now be slightly generalized; for example, a term 

£l • QbYI ■ h e 2 ■ Qb23 ■ £3 h • Gb31 • h 

would have to replaced by 

£l • Qb12 ■ k 2 S 2 - Gb23 ■ £3 h ■ Qb31 -fa - E\ ■ Qf\2 ■ k 2 £ 2 ■ 0F23 ' ^3 ^3 • Sf31 ' h 

Second, the above nonvanishing coincidence limits of Qb and Qp must be considered as cycles 
of length one, and included in the replacement rule fjl .4j) : 

QB{n,n) ->■ QB{n,n) - GF(n,Ti) (10.7) 

Third, there is now also a replacement rule for the determinant factor: 



det" 



sin(iJ) 



det" 



tan(Z) 



(10.8) 



This algorithm yields the one- loop N - photon amplitudes in a constant field for spinor QED. 

Although the worldline formalism in a constant field has already found extensive applications 
[52l E9 EH ESI [2U [561 EH EH [59] a systematic investigation of the IBP procedure for the constant 
field case has so far been lacking. Inspection of the symmetric partial integration algorithm of 
section d] shows that it applies as well as to this case, and it leads to essentially the same de- 
composition of the integrand into cycles and tails, with only two modifications: First, one-cycles 
have to be included; and second, the definition of the "bicycle of length n" , now to be denoted 
by Q{i\i 2 ■ ■ ■ i n ), has to be changed to 



Q(iii 2 ) := ^tr(Fjj • Qsi^ • F i2 ■ Qsnh) 

= £i • Qbi2 ■ k 2 e 2 ■ Qb2\ ■ k\ - e\ ■ Qb\2 • £2^2 ■ Qb2\ • h 
Q{i\i 2 ...i n ) := ^{F^ ■ g B hi 2 - F i2 ■ g B i 2 i A - ■ ■ F in - gBi n n) (n > 3) (10.9) 



These bicycles generalize the vacuum ones (|4.5|) . and will turn into them in the limit F — > 0. 

However, it turns out that the presence of one-cycles - which would lead to a significant prolif- 
eration of terms as compared to the vacuum case - can be altogether avoided, by the following 
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slight modification of the IBP: First, we observe that we have the freedom to shift both the 
Green's function Qbu and its derivatives Qbu by arbitrary constant matrices. For Qbu-i which 
appears only in the first term in the exponent of the master formula (110. ID . this is an obvious 
consequence of momentum conservation, ^ ■ kj = 0. For Qsiji we note that it can appear in a 
e i ■ Q Bi j - kj, ki-Q Bij • kj or £j • Q Bij • Sj ■ The first type of terms comes directly from the second term 
in the master formula, and here again a constant matrix added to Qbu drops out immediately 
because of momentum conservation. The second type of terms arises in the integration-by-parts 
procedure as a 



N N 

di ^2 -G B ijh - kj =}^ki- Qbu ■ kj 
l,j=i i=i 

with j running, so that again we can modify Qbu by a constant. The third type arises as the 
integral of a e% ■ Qbu • Sj in the IBP procedure, and here the constant matrix can be added as 
an integration constant. 

For the scalar QED case, we can use this freedom to directly eliminate all one-cycles by sub- 
tracting from Q B ij its (constant) coincidence limit, that is, replacing Q B ij by Qbu = QBij — Qbu 
throughout. For the spinor QED case, this would not make sense, since there are still the 
fermionic one-cycles. Here instead one should anticipate the application of the replacement rule 
for one-cycles (|10.8|) , and use the freedom of modifying Q Bij to replace it by 

Qbu '■= Qbu + Qfu = Qbu — Qbu + Gfu 
-izg B12 \ 

cot Z- tan Z\ (10.10) 

smZ J y ' 

This eliminates all one-cycles, since now by construction Qbu — Qfu = 0. 

With these modifications, we can write down the one- loop iV-photon amplitudes in a constant 
field for scalar and spinor QED in a way which is almost as compact as in vacuum. For example, 
the formula (|4.1ip for the four-photon amplitude in scalar QED generalizes to 



Tscal [kl , £\ ; • • • ; &4 , £4] 
r-T 



r (27r) D ,5(J> 



where 



(4tt) d / 2 

/ dndridrsdn (Q 4 4 + Q'l + 02 + Sf ) exp { V 



dT -m 2 T 1 i— A 

e det 2 



sin Z 



Q 4 4 

Q'i 



Q'i 2 



-ki ■ Qbu ■ kj 



0(1234) + 0(1243) + 0(1324) 

£7(123)7" (4) + 0(234)7" '(1) + 0(341)7" (2) + 0(412)7' (3) 
0(12)7'(34) + 0(13)7'(24) + 0(14)7"' (23) 
+0(23)7'(14) + 0(24)7'(13) + 0(34)7'(12) 
0(12)0(34) + 0(13)0(24) + 0(14)0(23) 



10.111 



(10.12) 
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Here the bicycles Q(i\ ■ ■ ■ i n ) differ from the original ones (|10.9p only by the replacement of Gb%j 
by Gsij, and the tails T{i\i 2 ■ ■ ■ i n ) are isomorphic to the one of the vacuum case (|4.6p . ([4.10p . 

a 

T'iij) ■= ^2 Si ■ QBia ■ k a Sj ■ GBjb ■ h + Ei ■ Gsij ' Bj (h • GlBia ' K ~ kj ■ Gsja ■ Kj 

a,b a 

(10.13) 

As before, a "prime" on a tail means that its cycle have been removed; note that now this is 
necessary already for the one-tail. 

The corresponding representation for spinor QED is obtained from ([10.110 . (|10.12p by the change 
of determinants (110. 8p . multiplication by the global factor of —2, and replacement of the bicycles 

G(h ... i n ) by the "super-bicycles" Gs(h ■ ■ ■ *'n) 

Gsihh) ■= (-^i ' ^Bhi 2 ■ F% 2 ■ GBi 2 i?j 

— -tr (F ix ■ GFhi 2 ■ Fi 2 ' GFish) 



Gsihh ■ ■ ■ in) '■— tr [F^ ■ GBni 2 • F i2 • GBi 2 i 3 • • • F in ■ GBi n ii 

— tr • GFhi 2 • F i2 ■ Gf% 2 v a ■ ■ ■ F in ■ QFinii) ( n - 3 ) 

(10.14) 

Moreover, the GBij's in the tails ([10.130 must also be replaced by Gsij's- 

The covariantization of the tails does not seem to extend to the constant field case in a natural 
manner, except for the modification of the one-tail ([3.80 . which generalizes to 



T rW :=2^ 1 (10.15) 

a 

Also the optimized form of the two-tail ([7.20 does not seem to generalize to the case of the 
general constant field. It does so, however, for the important special case of a self-dual field, 
which obeys F 2 ~ 1 [23 [581 [59] . Here one can generalize the total derivative (|7.ip to 



, ^^j2tr(FiF 2 )5i5 2 e ( ) + £\ ■ e 2 <9i<9 2 e ( ' ) - E\ ■ k 2 di (e 2 ■ Gb2J • fye^ - e 2 ■ hd 2 (ei ■ Gb 

and adding this to the (unsubtracted) two-tail T(12) of ([10.130 . one arrives at 

7ff(u) := ^J&a ■ GBai ■ Hij ■ GBjb ■ h (10.17) 



a,b 



with the same tensor as in ([7.30 
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11 Conclusions 



We have continued here the systematic investigation of the Bern-Kosower partial integration 
procedure, initiated in [36] and continued in [37]. We have presented an IBP algorithm that 
unambiguously leads to a form of the integrand of the one-loop N photon amplitude in scalar 
QED which is manifestly gauge invariant (transversal) at the integrand level, and suitable for 
an application of the Bern-Kosower rules. We have worked out this integrand explicitly at the 
four-point (two-tail) level. The N - point integrand contains N vectors ri,...,rjv which are 
constrained only by the condition (13.3H . Further study will be needed to find out what is the 
significance of this ambiguity, and how to make the best use of it. As far as concerns the on- 
shell N - photon/gluon amplitudes, one would surmise that the dependence of the integrand on 
the vectors T{ is related to the usual dependence on the reference vectors qi,... ,qN which one 
would normally have in the application of the spinor helicity formalism, but does not exist any 
more once all polarization vectors are absorbed into field strength tensors. And indeed, there 
is clearly a relation: for example, consider the case of the the on-shell A-photon amplitudes 
with all helicities positive. When using the P representation together with the standard spinor 
helicity formalism (see, e.g., |46] ) . one can remove all terms involving a Gsij^i'^j by choosing, in 
the spinor helicity formalism, the same reference vector qi = q for all legs, since then ef -e^ = 0. 
Similarly, in the S representation one could make disappear all the factors r% ■ Fi ■ Fj ■ rj by 
choosing all r j = r equal, and r on-shell, since then 



However, it is clear that this match between the freedom of choosing the r^s and the qi's cannot 
be a perfect one, since the r^s need not be chosen as on-shell. 

All our representations are valid off-shell. This makes them relevant for state-of-the-art calcula- 
tions already at the four-point level, since neither the four-photon nor the four-gluon amplitudes 
are presently available in the literature fully off-shell (for any spin in the loop). This fact is 
particularly conspicuous in the case of spinor QED, where the on-shell four-photon scattering 
amplitude was obtained already in 1951 by Karplus and Neumann [61] , and the extension to 
the case of two off-shell legs in 1971 by Costantini et al. [62J. Our integral representations for 
the QED four-photon amplitudes are, with any of the various definitions of the tails, manifestly 
finite term by term and thus suitable for a numerical evaluation as they stand. For analytical 
purposes one would still like to reduce the various parameter integrals appearing in them to 
scalar box, triangle and bubble integrals. This could be done using existing tensor reduction 
algorithms (see [63] and refs. therein), but in a companion paper [63] we will rather perform this 
tensor reduction in a way that is specifically adapted to the structure of the worldline integrals. 

Due to this validity off-shell our representations can also be used to construct, by sewing, all 
higher-loop N - photon amplitudes. From the calculation of the two-loop QED f3 function [T7] it 
is clear that when calculating those multiloop amplitudes in the worldline formalism significative 
simplifications can be expected from a judicious application of IBP. 





on account of the identity [60] 




(11.2) 
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The manifest transversality is probably a more significant issue in the nonabelian case. Here we 
have in mind not so much the calculation of on-shell gluon amplitudes, for which other extremely 
powerful methods have been developed in recent years (see, e.g., [65} 166] and refs. therein), but 
rather of the one-particle-irreducible off-shell N - gluon vertices, for which there is presently 
still a dearth of efficient methods. For the use of these amplitudes, e.g. in Schwinger-Dyson 
equations, it is usually important to have them in a form that separates them into transversal 
and non-transversal parts, which normally requires a tedious analysis of the nonabelian Ward 
identities. This is one of the reasons why, using standard methods, the explicit calculation of 
these off-shell vertices has been completed so far only for the N = 3 case [391 167L f5T) . In another 
companion paper [35] we use the S representation to recalculate the three-gluon vertex, for the 
scalar, spinor, and gluon loops, achieving a drastic reduction in computational effort compared 
to earlier attempts (a summary of this calculation was given in [68J). The main advantages 
of our approach are that the gluon and spinor loop cases can be effortlessly obtained from the 
scalar loop one through the (off-shell extended) Bern-Kosower replacement rules, and that there 
is no necessity to solve the Ward identities, rather the decomposition into transversal and non- 
transversal pieces emerges automatically in the IBP procedure. In relation with the latter point 
it must also be mentioned that the boundary terms in the IBP procedure applied to the N- 
vertex always involve color-commutators, and are always connected to some lower-point term, 
even appearing with the same integral. Thus the algorithm makes it also easy to separate the 
genuinely new structures appearing in the N- vertex from those that, in terms of the effective 
action, only serve the completion of lower-point expressions to fully gauge-invariant ones; see 
05]. 

In the abelian case, we have also extended the IBP procedure to the QED iV-photon amplitudes 
in a constant external field. Here our main motivation is that these amplitudes can be used 
for the construction of higher-loop Euler-Heisenberg Lagrangians, and those as a tool to obtain 
insight about the asymptotic behavior of the QED N - photon amplitudes at high loop orders 
and large photon number [69\ [58] [59] [70] [71] . In that context it would be highly desirable to 
calculate the three-loop Euler-Heisenberg Lagrangian in various dimensions, and indeed the IBP 
procedure presented in section [10] has made it possible to finally achieve this goal at least for 
the 1+1 dimensional case [72] , 

It should be possible, and very interesting, to generalize our approach to the inclusion of grav- 
ity. On-shell, the gauge theory Bern-Kosower rules were generalized to the construction of 
string-based representations for the one-loop iV-graviton amplitudes in [73], and these gravity 
rules were then successfully applied at the four-point level in [73]. They also involve an IBP 
and "replacement rules" connecting the amplitudes with different spins in the loop. Worldline 
path integral representations of the one-loop effective actions for gravity have so far been con- 
structed for spin zero [24], spin half [25] [26] and spin one [27] in the loop, and can be used to 
obtain parameter integral representations of the corresponding off-shell one-particle irreducible 
iV-graviton amplitudes that are closely related to those string-based representations. They are 
written in terms of the same worldline Green's functions as the ones for gauge theory, and the 
challenge is again to find an IBP algorithm that would allow one to apply the replacement rules 
and at the same time make covariance manifest, where the latter now means the emergence of 
full Riemann tensors in the IBP. This algorithm can, however, not be a simple extension of the 
one which we have presented here, as one can see from the fact that in gravity there are no 
boundary terms in the IBP, so that the nonlinear terms in the Riemann tensor now have to be 
created by 5- functions; therefore a complete removal of all G_b's in the IBP is not called for, and 
in fact also not possible, as one can easily see already from the case of the graviton propagator. 



23 



Finally, in open string theory the one-loop gauge boson amplitudes can be written in terms 
of a master formula that is analogous to the master formula (jl.ip . only that the variables Tj 
parametrize positions along the boundary of the string worldsheet, and the Green's functions are 
worldsheet Green's functions (see, e.g., [75]). Since in all of our manipulations we have, apart 
from the translation invariance in proper-time, not used any specific properties of the worldline 
Green's functions, our IBP procedure could as well be applied at the string level to achieve a 
form factor decomposition based on gauge invariance. 
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